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Abstract. The variational principle and the corresponding differential equation for geodesic 
circles in two dimensional (pseudo)-Riemannian space are being discovered. The relation- 
ship with the physical notion of uniformly accelerated relativistic particle is emphasized. 
The known form of spin-curvature interaction emerges due to the presence of second or- 
der derivatives in the expression for the Lagrange function. The variational equation itself 
reduces to the unique invariant variational equation of constant Frenct curvature in two 
dimensional (pseudo)-Euclidcan geometry. 
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1 Introduction 

It turns out that the notion of the uniformly accelerated relativistic test particle world line [I] 
coincides with the notion of the geodesic circle in pseudo-Riemannian geometry [2]. In two 
dimensions both these notions may be equivalently replaced by the only condition that the first 
(and the only) Frenet curvature k of the curve in the consideration keeps constant along this 
curve (in physical terms, along the world line of the test particle). 

In the natural parametrization by the arc length s (in physical terms, by proper time) all 
the notions, mentioned above, amount to the following third order differential equation (cf. 
Appendix lA.Sp : 

jjS^l j-)2^m jj2^n j^^J ^ 

ds^ ds"^ ds"^ ds 

The left hand side of this equation is known in the physical literature under the name of the 
Abraham vector. This vector is believed to adequately represent the notion of the relativistic 
acceleration of test particle [T] . Of course, the notion of geodesic circles does not depend on any 
arbitrary reparametrization s(.^) of the independent variable along the curve in (II. ip . 

One may ask whether geodesic circles could be the extremals of some variational problem. 
Since the equation (jl.ip involves third order derivatives, it is appropriate to speak of higher order 
variational calculus. Rather trying in naive way to find a Lagrange function for the equation (II. ip 
(or, better, an equivalent equation), it may happen more efficient to set a general problem of 
finding all possible third order variational equations in two-dimensional space which keep the 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
ihttp: / /www. emis.de/journals / SIGMA /symmetry2007.html | 
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Frenet curvature k constant along their solutions. Of course, one may hope to solve in reasonable 
fullness such an inverse problem only if restricted to certain types of the space geometry. The 
natural restriction with regard to (pseudo)-Riemannian space should, of course, be that of its 
local model - the (pseudo)-Euclidean space. Thus one comes to the necessity of the formulation 
of the invariant inverse variational problem possessing the first integral k. This problem in 
two dimensions received the ultimate answer in [3] which is quoted below as Proposition [2l 
This proposition says that to obtain geodesic circles, one should, in two dimensions, necessarily 
handle the Frenet curvature itself as the integrand of the variational functional. 

The fashion of including Frenet curvatures into the variational integrand in physics literature, 
initiated, in our opinion, by the papers of Plyushchay and Nesterenko, counts 30 years by 
now. Still, sometimes it escapes the common knowledge that such quantities as, say, Frenet 
curvatures, are always constants of motion along the extremals of the variational problems where 
the Lagrange functions depend on nothing but the corresponding curvatures themselves (4]. This 
fact being noticed long ago [5] for the (pseudo)-Euclidean case, we present it here in general 
framework by means of rather trivial observation of Proposition [TJ Unfortunately, only in 
two dimensions the first Frenet curvature taken as the Lagrange function produces, due to its 
linearity in second derivatives, a third order equation. 

In Section 13.21 we deduce the variational equation from the variational functional 

J kdi + ds 

in (pseudo)-Riemannian two-dimensional space and prove that the set of its extremals includes 
all geodesic circles. 

Also we discovered quite an interesting fact, in our opinion, that the force on the particle, 
having the same form as that produced by spin-curvature interaction, emerges as a simple result 
of calculating the variation of second derivatives in the variational integrand. 



2 Preliminaries 

2.1 Parameter independence and Ostrohrads'kyj mechanics 

The right action of the prolonged group GL(^2)(^) =^ <^^o(I^)I^)o of parameter transformations 
(invertible transformations of the independent variable ^) on the second order velocities space 
T^M = {x^jU'",^"} gives rise to the so-called fundamental fields on T^M (cf. Appendix lA.ip : 

A function / defined on T^M does not depend on the change of independent variable 
(so-called parameter-independence) if and only if 

Ci/ = 0, C2/ = 0. (2.2) 

On the other hand, a function L defined on T^M constitutes a parameter-independent vari- 
ational problem with the functional J L{x'^,u",u")dS^ if and only if the following Zermelo con- 
ditions are satisfied [U formula (8.19)]: 

CiL = L, C2L = 0. (2.3) 



Note 1. At this point it worths mentioning the difference between the conditions p. 21) and ()2.3p . 
It would be instructive not to confuse the idea of the parameter independence of a variational 
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problem with that of the parameter independence of Lagrange function. The Lagrange function 
of a parameter-independent variational problem in not parameter-independent itself. 

The formulation of Zermelo conditions in the most general case of arbitrary order and of 
arbitrary number of independent variables (Zermelo-Geheniau) may be read from [7j whereas 
the corresponding proof may be found in [S] as well as in [3]. 

Let us introduce the generalized momenta: 

The Euler-Poisson equation is given by 

£n = ^-^ = 0. (2.5) 

The Hamilton function is given by: 

H = p^^^u'' +PnU'' -L. 
Lemma 1 ([3 [9]). 

di 

Proposition 1. If a function Ln is parameter-independent and if a function Li constitutes a pa- 
rameter-independent variational problem, then Lji is constant along the extremals of L = Ln+Li. 

Proof. By Lemma [1] and in force of the properties (j2.2p and (j2.3p we calculate i/^jj+^j = 
Ci(Ln + Li) — -^C2{Lii + Li) — L = —Lu. But as far as the Hamilton function is constant of 
motion, so is the Ln- ^ 

2.2 Inverse variational problem for geodesic circles 

As stated in the previous papers by the author, the general third-order variational equation 
takes on the shape [TT] : 

AmnU'^ + U^d^i AmnU"" + ^mn^" + Cm = 0, (2.6) 

where the skew-symmetric matrix A, the matrix B, and a column c all depend on x", and n", 
and satisfy the following system of partial differential equations: 

2 5 - 4 (9 r^^n] + (9 + 2 9 fc = 0, 

2 d^k d^lmCn] - 4 l9^[m^n] k + '^i'^ dj. Amn + 6 dj^m^nk] = 0, 
4<9^[mC„] - 2D^(9^[^C„] - Y)^^ Amn = 0. 

Here the differential operator is the lowest order truncated operator of total derivative, 
Dj^ = d^ + vl'd^i. We refer to App endix I A . 2 1 on some additional comments concerning the nature 
of the variationality conditions (|2.7p . 
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On the other hand, we are interested only in differential equations which enjoy the (pseudo)- 
Euclidean symmetry. Because of that in two dimensions the skew-symmetric matrix A in (|2.6|) 
is invertible, there is no problem to formulate the concept of the symmetry in the way 

{X£n)\s^^ = 0, (2.8) 

where £ denotes the left hand side of the equation (j2.6p and X stands for the infinitesimal 
generator of the (pseudo)-Euclidean group. 

The system of partial differential equations (j2.7p together with (j2.8p was solved in [3] and 
the following result was established: 

Proposition 2. Let some system of third order differential equations ()2.6p enjoy the following 
properties: 

• the conditions ()2.7p are satisfied; 

• the system (12. 6p possesses Euclidean symmetry according to ()2.8p ." 

• the Euclidean geodesies u = enter in the set of solutions of (|2.6p ; 
dk 

• — = along the solutions of (12. 6p . 
Then 

e„iil^ (u-u) J (u-u)un — (u-u)un , s 

^" = T^-3V^f"'^ +"^^ ^^^^ —■ 2.9 

In the above statement the denotation e„m stands for the skew-symmetric Levi-Civita symbol. 
This expression (12. 9p may be obtained as an Euler-Poisson expression for the Lagrange function 

r Cmn'W ^ ii ii /n in\ 

L = 3 m \\u\\. (2.10) 

11^11 

The first addend in (j2.10p is sometimes called the signed Frenet curvature [12] in E^. This, along 
with the observation that in two dimensional (pseudo)-Riemannian space the Frenet curvature 

r^ = ±v^ " „3 (2-11) 



WuW \\u\ 



depends linearly on the covariant derivative u' and thus produces at most the third order Euler- 
Poisson equation, suggests the next assertion, based on Proposition [TJ 

Proposition 3 (|131 [3]). The variational functional with the Lagrange function 

tTZ /n ffcn'"'^'"'" II II ir,\ 

L = \/\g\ — |j — m \\u\\ (2.12) 

produces geodesic circles in two dimensional Riemannian space. 

Note 2. The Lagrange function (j2.10p looks like a relativistic analogue of the one recently 
treated in 
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3 The variational equation for geodesic circles 
3.1 The generalized covariant momenta 

It is known that the Euler-Poisson expression (in this paper, - of the third order) (j2.5p con- 
stitutes a covariant geometric object. But the conventional momentum pn from (j2.4p does 
not. Therefore we introduce covariant momenta tt and tt^^^ that, in the case of the Lagrange 
function we shall deal with, represent some relative vectors in the following manner. First, we 

assume that the Lagrange function depends on the variables it" exclusively through the covariant 

dL dL 

derivative u'. Moreover, in the case handled in this article, the partial derivatives — — and — 

9m" 5u'" 

coincide. Let us introduce the "truncated" partial derivatives 
dL dL du'i dL 

(3.1) 



Ou" du^ du^ du't 
and, for future use, 

dL _ dL du'i dL 
Then the quantities 

= 1^ and (3.3) 
du 

dL 

TT = — - 7r(i)' (see ^) (3.4) 
du 

flT 

= — - 2 - 7r(i) (by virtue of ^) (3.5) 

may each be applied the operation of the covariant differentiation according to the formu- 
lae ()A.5p . For example, the covariant momentum vr^^^ produces its covariant derivative 

^(')'n = ^(')n-r™znvr«„n'. (3.6) 

In the above terms the quantity pn of (j2.4p is expressed as follows: 

= 7r„ + 2 V^nu'^TT^^^ + ^^^^ 'n (on basc of ([33])) 

- -fr^^ln (on base of (ISTSD ) 

= 7r„ + rVnu'^vrW (on base of (IMI)). (3.7) 

Differentiating (j3.7p and applying (jA.SP to express the ordinary derivatives of 7r„ and in 
terms of tt' and u' along with (j3.6p . gives: 

f)-pi 

Pn = (^'n + r'™„vr,n-) + ^L^n'u^^TrW, 

+ (r'^„n'- - r'^„r-,fcn%'=)7r(i), + r'^„n™(7r(i)', + ThiTT^^^u") 

dT% 

dx^ 

Now the Euler-Poisson expression ()2.5p may be handled with the help of our "truncated" x-de- 
rivative (j3.2p in the following way: 

dL 

£n=Q^-^'n- {tT^^'^'i + VrOr'™„n- - 7r('\T' „,nu' - TT^'^U^ u'' R^kJ , (3.8) 

making use of the definition of the curvature tensor ()A.9p . 



+ 7r^^on U \Lmn^^lk-\ TT-r. i-^lni^ mk 
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3.2 The Euler— Poisson expression for the signed Prenet curvature 

Let us carry out the previous considerations taking for the Lagrange function Ln the signed 
Frenet curvature: 



\9\ 



\U\ 



(3. 



Making use of (|A.7|) and (jA.Gh one gets immediately 



W\9\ 



The expression for 7r(^) is obtained from ()3.3p in still easier way: 

,k 



(1) 



u 



(3.10) 



(3.11) 



The covariant derivative of 7r(i) is presented by the decomposition formula (j3.6p . While calcu- 
lating ^(^)„, in (|3.6p we profit both from ()A.7p and from (jA.Sp . Then we replace 7r(^)„, in (|3.6I 
with (13. lip to obtain: 



vr 



(1)' 



{u-u')eknvf' 



\9\ 



k„,m 



\\U\\ \\U\\ 

Finally, the simplification formula ()A.4p may be used to produce 



n 



Ik 



{u-u')eknu'' 



\\u\\ \\u\\ 

The calculation of the sum vr + tt*-^)' on grounds of (j3.4p is straightforward: 

,lk 



(3.12) 



(3.13) 



Ifl 



UneikU^u"' 



\\u\\ \\u\ 
Or, alternatively, utilizing (IA.3p . one obtains 



(1)' 



W\9\ 



[u-u') eknu'' 



\\u\\ \\u\ 
To calculate tt, we extract (j3.13p from p.l4p : 



151 



/fc 



It 



(3.14) 



(3.15) 



Now, the Euler-Poisson expression (|3.8p for the Lagrange function (|3.9p by virtue of (|3.10p . 
(f3TT]l . and (IXTil) reduces to 



^n(-^^Il) = -T^'n - T^^^^Rn. 



(3.16) 



where the covariant derivative of the momentum tt may be calculated starting with the defini- 
ton (jA.SP applied to (j3.15p . wherein (u'^)' should be calculated again along the lines of (jA.SP : 



/k\ 



U 



Ilk 
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Using then ()A.7P and (jA.SP along with the simplification formula (|A.4p once more, we calculate 

= -<5 ;--7^ ^knU + - — CknU ■ 

\\u\\ 

Thus in the disclosed form the Euler-Poisson expression for the total Lagrange function (j2.12p 
reads: 

'^n — V \9\ — T3 roy'\g\ - -g 

\\u\\ \\u\\ 

{U-U)u'n- {u' •u)Un \/M q m I k /o 1 

+ m ——n -^—n enigRnki^u" u u . (3.17) 

||it|| 

The Euler-Poisson equation, solved with respect to highest derivative, now is 



IX w / \\u\ 



where we used the notion of the contravariant skew-symmetric Levi-Civita symbol e^' 
<^^^[9pq\ 9^^ (^kh ™d introduced the shortcut 



^« = ^mgRnkl''u"'u^u'' (3.19) 

\\u\\ 

to denote the force on the particle, evoked by the curvature of the (pseudo)-Riemannian struc- 
ture. Equation ()3.17p obviously generalizes that of (12. 9p . 



Note 3. If, on the other hand, we had introduced a spin tensor Smq = T^^^^mUq — T^^^^Um, 
the second term on the right of formula (|3.16p would have got an interpretation as the force, 
evoked by the existence of quasi-classical spin, and would coincide with that present in Dixon's 
equation [15]. 



Proposition 4. The differential equation (I3.18P describes all geodesic circles. 

Proof. We have to show that, for each solution of (jl.ip it is always possible to perform such 
transformation of the independent variable s to ^ that after that change of variables the equa- 
tion ([338]) will hold. 

First, let us rewrite the differential equation of geodesic circles (jl.ip in arbitrary paramet- 
rization: 

u" u-u" iu-u') , iu-u'f 
5" "" + 3^ 5^" Ji-u. (3.20) 



We add to it one more equation, which will play the role of such a one, that fixes the parametriza- 
tion along the curve: 

u-u" {u-u'Y\ „h f m , ^ \ , 

This equation is consistent with the equation ()3.18p because it presents nothing more than a 

mere consequence of (|3.18p . obtained by means of the contraction with u^- Further, let us 
yii k 

substitute o in ()3.18p from the equation (I3.20p : 

llwll 
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It remains to insert the left hand side of (j3.2ip into the left hand side of (|3.22p and to notice 
that in two-dimensional case what comes out is the identity: 

J_ fjILe^'u'nUi - e^^ui-JZn) = e'=" f ^ ^'^ - m ^ n„ - . (3.23) 
u-u \\\u\\ J \\u\\ J 

To see that (|3.23p is satisfied identically, it suffices to apply the simplification formula ()A.3P 
separately to the terms involving IZ and then to what remains, and to notice that u^TZn = 
due to the previously introduced notation (|3.19p . ■ 



A Appendix 

A.l Parameter independence and fundamental fields 

The T^'M: loose comments. In the variational calculus various objects depend not only on 
the variables of the configuration space, but also on the velocities and on the accelerations of 
the first as well as of higher orders. The space of this extended number of variables may be 
introduced in different ways, as it has in fact been by different authors. In present paper we 
choose to use the definition, belonging to Ehresmann [T7], of the higher order velocity space 
T''M = J''o(K, M)o as the set of jets of mappings from a neighborhood of the origin in M to the 
configuration space M = rr", which all start at the origin G M. This space T^M is naturally 
endowed by several geometric structures. The Reader may consult the monographs [H] and [9] 
on the subject. For the purposes of this paper it suffices to think about T'^M as a manifold 
of the variables x*^, u", . . . , Ur-ii constructed of the successive derivatives of the configuration 
space variables by the independent evolutionary one. 

The fundamental fields. The group of invertible jets J^o(I^)I^)o with the source as well 
as the target at the origin G M presents an appropriate geometrical concept when speaking 
about local transformations of the parameter along the curve in a manifold. Let again the 
variables = x", = x" denote the first and the second derivatives of coordinates along the 
curve x"(^), so that the jet Jq x at zero is presented by the array of second-order polynomials 
it"^ + n"^^. Another jet, j'-^-'o" S J^o(1^5 1^)01 which is presented by the polynomial + 
acts on the right upon the previous one by the composition: 

= vJ'ai + \{u''P + u''a^)f modo(e^). (A.l) 
Consider an e-shift (JeiC) of the transformation a of the independent variable ^ in the local curve 
expression x"(^). This shift is presented in J^o(lR)IR)o by some = — ^(0) and (5e = — ^(0). 

It evokes a corresponding flow {u", it^} € T^A/ induced by the action (jA.ip : 

The generator of this flow is: 

+ ( '^yn + 2a,^u 



de 



ou"- \ de de 



e=0 



' (A.2) 



=0 

Now it suffices to limit ourselves to the shift cr£(^) of the form 0"e(^) = £, + eT{(,). Then the 
generator (jA.2p reads: 

„ dr , . d / „ d'^T , . „ dr , .\ d 
« ^ 0) 7^ + U" —^ 0) + 2'it" —(0 — . 

Taking r(^) = ^ we obtain the first fundamental field d in (j2.ip . Choosing t(^) = we 
obtain the second fundamental field C2 in (|2.ip . 
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A. 2 The inverse variational problem: loose comments 

The system of equations ()2.7p along with the specific guise ()2.6p of the Euler-Poisson equation 
arises as the general solution of the so-called inverse variational problem in the calculus of 
variations. This problem has been attacked by numerous authors from varying points of view 
based on different approaches. Roughly speaking, the problem consists in finding out the 
criterion that for an a priori given differential equation there locally exists a Lagrange function, 
from which this equation follows by the variational procedure, applied to the corresponding 
action functional. As far as the act of the variation of the functional may be expressed in the 
form of some operator 6 action on the Lagrange function L in such a way that £ = 5L he 
a well-defined geometric object, represented in local coordinates by means of well-known system 
of the Euler-Poisson expressions it is tempting to give such a definition of 6, wherewith 
the cohomology complex property 5^ should hold. The equations (12. 7p express the fact that the 
differential form E^dx^ is closed with respect to 5. It is not our intention in this paper to discuss 
further the ways of defining the operator 5. Along with the sources |10l [TT] we wish to show the 
interested reader to the book by Olga Krupkova [18] with the plentitude of references therein. 
However, our explanation here was based on the approach of Tulczyjew |19j . 



A. 3 Simplifications of exterior products 

In two dimensions some vector and tensor skew-symmetric expressions simplify drastically. 
Let a", fe" and c" denote arbitrary vectors and let T^mn denote for a moment an arbitrary 
three-index quantity. The following two simplifications keep true if the underlying (pseudo)- 
Riemannian manifold is two-dimensional: 

gm.nO^a'^eikb^ - granO^b'' eikO^ + OkemnO^b'' = 0, (A.3) 
^mnCl I ikC - emnb I IkO C + Cmna 1 IkO C — U. lA.4j 

The proof consists in the ingenuous calculation. 



A. 4 Formulae from (pseudo)-Riemannian geometry 

Here we list some well-known relations, involving Christoffel symbols, curvature tensor, and the 
covariant derivative. The latter will be denoted by the prime superscript, while to denote the 
ordinary derivative the dot will be used, e.g. n" = x", it"" = x", etc. 

a' - = + T^ima'^u', a'n = d„ - r™/„a^u', (A.5) 

Q^k ~ Sml^ kn + 9nl>- km, V^-OJ 

bl = vM r'in, where g = dei[gnm], (A. 7) 



dx"' 
( 1 



\u\\ / It 



|5 ■ 



(A. 



„ t _ fcn t/1 rnn ™ _W -pg /. 

^'-kmn — Qx"^ dx^ "i<? kn kq'- mn- V-*^' 



A.5 Characterization of geodesic circles 

It is not difficult to check that the condition ^ = in two dimensions is equal to the equa- 
tion (jl.ip on the shell gmnU^u^ = 1- More precisely, at the constraint manifold 



u - u = 1, 



(A.IO) 
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uu' = 0, (A.ll) 
u' ■ u' + u ■ u" = 0, (A.12) 

the differential equation 

k^' = ±u'-u" = (A.13) 
is equivalent to the differential equation of the geodesic circles 

u" +{u'-u')u = 0. (A.14) 
In fact, from ()A.13P and ()A.12p one solves for u": 

U I = ; U -U , 

and then implementing ()A.3P with the help of (|A.10P and (|A.lip obtains ()A.14p . Conversely, 
the (jA.lSp in nothing but the (|A.14p contracted with u' at the constraint manifold (|A.lip . 
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